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FOR  EWORD 

A systematic  study  of  spin  - stabilized  rockets  at  the  U.  S. 
Naval  Ordnance  Test  Station  has  provided,  among  other  things, 
pitching  moment  and  normal  force  coefficient  data  on  short  bodies 
of  revolution.  These  data  were  for  5 , 7-  and  9-caliber  body 
lengths  at  a Mach  number  of  0.  26,  and  included  runs  with  and 
without  spin.  An  attempt  to  explain  the  nonlinear  coefficients 
fnr  angles  of  attack  up  to  30  degrees  and  with  zero  spin  has 
resulted  in  the  method  presented  here.  The  success  of  this 
theory  at  low  speed  prompted  its  extension  to  high  speed  flow. 

This  study  has  been  supported  by  NOTS  Local  Project  LP  701. 

This  report  has  been  reviewed  for  technical  adequacy  by  Dr. 
Henry  T.  Nagamatsu,  personal  services  contractor  to  NOTS,  of 
the  California  Institute  of  Technology,  Pasadena. 

4$. \ 
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He^, /Aerodynamic s Branch 

This  memorandum  is  transmitted  for  information  purposes 
only.  It  does  not  represent  the  official  views  or  final 
judgment  of  the  Naval  Ordnance  Test  Station,  and  the 
Station  assumes  no  responsibility  for  action  taken  on  the 
basis  of  its  contents. 
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ABSTRACT 

A method  in  developed  which  accurately  predicts  for  blunt- 
based  bodies  of  revolution  the  normal  force  coefficient  and  the 
pitching  moment  coefficient  for  angles  of  attack  far  beyond  the 
range  of  potential  theory.  This  method  is  based  on  the  principle 
of  superposition  of  the  results  of  potential  theory  and  the  viscous 
force  on  a cylindrical  body  due  to  the  transverse  component  of 
flow.  In  contrast  to  previously  used  methods,  the  viscous  cross 
force  is  assumed  not  to  be  in  a steady  state,  but  in  a transient 
development  along  the  body. 

This  method  is  compared  with  experimental  data  for  both 
subsonic  and  supersonic  flows  and  with  both  laminar  and  turbulent 
axial  boundary  layer.  The  deviations  between  theory  and  experi- 
ment are  in  most  cases  an  order  of  magnitude  om^ller  than  for 
a previous  method  due  to  Allen  and  Perkins. 
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INTRODUCTION 

The  aerodynamic  design  of  aircraft  rockets  is  somewhat 
unique  in  that  the  most  practical  body  shape  from  other  design 
considerations  usually  is  a cylindrical  shape  with  a blunt  base 
and  a moderately  blunt  nose.  Very  little  good  theoretical  infor- 
mation is  available  on  the  lift  and  pitching  moment  for  bodies  of 
this  type,  except  at  very  small  pitching  angles.  For  this  reason, 
a considerable  number  of  wind-tunnel  tests  have  been  made  by 
the  Naval  Ordnance  Test  Station  to  build  up  a fund  of  information 
for  design  purposes. 

While  examining  some  of  these  data  for  the  purpose  of 
correlation  between  tests  and  with  available  theory,  it  was  found 
that  certain  subsonic  test  data  on  short  bodies  at  large  angles  of 
attack  differed  widely  from  any  theoretical  predictions  available. 
The  theory  of  Max  Munk  (Ref.  1)  waa  found  to  give  good  results 
at  small  angles;  this  is  satisfactory,  since  his  is  a linearized 
potential  theory,  and  is  not  expected  to  be  valid  for  large  angles. 
A theory  proposed  by  H.  Julian  Allen  (Ref.  2),  which  has  enjoyed 
a measure  of  success  for  supersonic  flow  data,  grossly  over- 
estimates the  effects  of  viscous  crossflow  at  low  subsonic  speeds, 
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where  it  would  at  first  be  expected  to  be  at  its  best. 

The  theory  of  Allen  has  been  carefully  examined,  and  a 
completely  new  version  of  his  method  is  presented  here.  The 
problem  of  viscous  crossflow  has  been  reformulated  with  results 
quite  different  from  his.  His  universal  use  of  Munk's  results  for 
the  linear  potential  theory  contribution  is  questioned,  and  a 
revision  of  this  policy  is  recommended.  The  results  of  a number 
of  wind-tunnel  tests  by  NOTS  are  presented  in  support  of  this  new 
theory,  showing  excellent  correlation.  The  predictions  of  Allen's 
theory  are  included  for  comparison,  and  the  results  of  Munk's 
linear  theory  included  as  a sort  of  reference  line,  without  in- 
tending a comparison. 

This  method  is  found  to  give  formulas  for  certain  cases  that 
would  appear  to  make  it  more  difficult  to  use  for  engineering 
design  purposes  than  Allen's  method.  It  is  shown,  however, 
that  considerable  simplification  can  often  be  effected  for  this  type 
of  application,  without  seriously  impairing  the  accuracy.  An 
important  result  of  this  method  is  that  it  provides  a means  of 
extrapolating  the  results  of  tests  at  low  Reynolds  numbers,  with 
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completely  laminar  boundary  layer,  to  high  Reynolds  numbers, 
with  fully  developed  turbulent  boundary  layers.  This  is  a 
valuable  asset  for  design  work  when  test  facilities  are  limited. 

THEORY 

Historical  Development 

The  classical  work  of  Munk  (Ref.  1)  is  now  so  well  known  that 
the  results  of  his  theory  may  be  quoted  without  detailed  explana 
tion.  He  found  that  the  lift  force  distribution  along  the  body  of 
an  airship  is  given  by 

dF  = (k2  - k,  )~  sin  2adx  (l) 

which  integrates  to  zero  along  the  body  of  the  airship,  but  for  a 
body  of  revolution  with  a flat  base  of  area  Sg,  gives 

L = J*dF  = (k2  - k,  )Sfi  sin  2a  (2) 

For  a long  slender  body,  k2  k,  is  almost  unity,  and  if  there  is 
no  boattail-'So  that  the  base  area  Sg  is  the  reference  area  for 
defining  lift  coefficient- -this  coefficient  is  usually  approximated 


as 
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CL  = 2a 


(3) 


Munk  realized  that  his  linearized  potential  theory  would  apply 
only  at  small  angles  of  attack,  and  that  it  would  fail  at  large 
angles,  chiefly  due  to  viscous  effects.  There  was  apparently 
little  interest  at  that  time  in  developing  a theory  for  large  angles. 
Recently  Allen  (Ref.  2,  3,  4)  has  developed  a method  for  pre- 
dicting lift,  pitching-moment,  nd  the  drag  increment  due  to  yaw, 
that  has  enjoyed  considerable  success  with  slender  bodies  of 
revolution  at  supersonic  speeds. 

Allen  has  assumed  that  the  viscous  contribution  to  lift  is 
independent  of  the  contribution  of  the  potential  flow,  and  writes 
the  lift  coefficient  and  pitching  moment  coefficient  as  the  sum 
of  two  terms : 


SB  A, 


Cl  = (k  j - k,  )-^-sin  2 a cos-^-  + i)Cp  -j-sin2acosa  (4a) 


CM  = t 


V - SR(i  - xm)  , 

] ( k 2 - k,)sin  2acos-2_ 


+ nCDc*X-  (xm  - xc)sin2a 


(4b) 
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where  the  first  term  is  always  the  result  of  Munk's  theory, 
modifted  by  the  results  of  Ward  (Ref.  5),  v/ho  showed  that  the 
potential  cross  force  is  directed  midway  between  the  normal  to 
the  axis  of  revolution  and  the  normal  to  the  wind  direction. 

Drag  coefficient  increments  are  not  being  treated  in  this  report, 

■’o  Allen's  drag  formula  is  omitted.  The  quantity  Cp  as  used  by 
Allen  refers  to  the  steady-state  drag  coefficient  of  an  infinite 
cylinder  in  a flow  normal  to  its  axis.  This  was  found  by  Allen 
to  be  Cj)c  ; 1.2  for  a wide  range  of  Reynolds  numbers  correspond- 
ing to  the  usual  range  of  crossflow  Reynolds  numbers  being  con- 
sidered. The  factor  T]  corrects  this  cylinder  crossflow  drag 
for  the  finite  finenet  s ratio.  Values  of  q have  been  tabulated 
by  Goldstein  (Ref.  6).  The  xm  is  the  (arbitrary)  point  about 
which  pitching  moments  are  measured;  xc  is  the  axial  position 
of  the  centroid  of  the  planform  area  Ap,  and  the  factor  sin  a 
arises  from  the  transformation  of  dynamic  pressures  from  the 
cross  velocity  V sin  3 to  the  free  Btream  velocity  V . 
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Since  his  theory  is  only  approximate,  Allen  simplified  his 
formulas  by  making  the  approximations  sin  <3  = a,  k2  k,  = 1, 
cos  a = 1,  This  results  in 


SB  Ar 

CL:2(y)at  nCDc(-^)«‘ 


(5a) 


CM  = *1' 


v - Sb(-^-  *m> 


la  ♦ 1CDc(-jf-)(Km  *c)aJ  (5b) 


and  his  theory  is  used  in  essentially  this  form,  J.  A.  F.  Hill  at 
the  Naval  Supersonic  Laboratory,  Massachusetts  Institute  of 
Technology,  used  Allen's  results  in  1950  (Ref.  7),  but  modified 
them  to  include  boundary  layer  displacement  thickness  in  the 
potential  theory  term,  thus  changing  the  effective  values  of 
both  Sg  and  V . He  also  observed  that  Allen's  assumed  value 
of  = 1.2  is  not  appropriate  if  the  axial  boundary  layer  is 

turbulent,  even  though  the  crossflow  Reynolds  number  corresponds 
to  a laminar  boundary  layer. 
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1 he  Viscous  Cross  Flow 

In  the  theory  presented  in  this  paper  we  shall  replace  the 
use  of  the  lift  coefficient  C^,  by  the  normal  force  coefficient 
Cpj  This  will  not  change  the  first  (potential  theory)  term  of 
equation  4a,  and  will  eliminate  the  factor  of  cos  a in  the  second 
(viscous)  term  of  this  equation.  It  also  eliminates  the  necessity 
of  Allen's  assumption  that  the  component  of  axial  drag  involved 
in  the  lift  force  may  be  neglected.  This  assumption  would  be 
invalid  for  the  angular  range  of  interest  here.  As  Hill  has  pointed 
out  in  Ref.  7,  Allen's  equation  (5a)  would  be  more  appropriate 
for  the  normal  force  coefficient  Cjq  for  comparison  purposes  here. 

Allen's  fundamental  equation  for  the  viscous  cross  flow  is 
apparently  quite  correct  if  properly  interpreted.  He  writes  the 
viscous  contribution  to  the  cross  force  from  a cylindrical  element 
of  length  dx  as 


dF  = 2rCD 


2 


dx 


(6) 
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where  r is  the  body  radius  at  the  point  x,  and  Vc  is  the  cross 
velocity,  as  shown  in  Fig.  1.  The  quantity  p is  the  usual  mass 
density,  and  Cj->c  the  drag  coefficient  of  a circular  cylinder  at 
the  Reynolds  number 

2rVc 

Rec  = — (7a) 


and  the  Mach  number 


(7b) 


with  the  kinematic  viscosity  v and  speed  of  sound  a evaluated  in 
the  undisturbed  flow.  Then,  since 


VL  = U «in  a 


it  follows  that 


a 


(7c) 


(8) 


in  terms  of  the  dynamic  pressure 


of  the  free  stream. 
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Since  the  drag  coeificient  of  a circular  cylinder  in  a cross 
flow  is  usually  given  for  the  two-dimensional  case  of  an  infinite 
cylinder,  Allen  adds  a parameter  , which  represents  the 
reduction  in  drag  coefficient  of  a cylinder  due  to  its  finite  length. 
Then  if  the  moment  coefficient  for  a body  of  diameter  d,  referred 
to  a point  xm, 


(9a) 


and  the  normal  force  coefficient 


C 


(9b) 


we  find  the  increments  in  these  coefficients,  due  to  viscosity, 
to  be 


A r _ 2 H s»n' 

acn — 


dx 


(10a) 


and 


AC 


M 


x)dx 


(10b) 
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which  lead  to  the  final  terms  of  equations  (4a)  and  (4b).  These 

terms  are  entirely  due  to  viscous  effects,  since  Cn  always 

c 

vanishes  for  a true  potential  flow. 

Evaluation  of  Parameters 

The  departure  of  the  theory  of  this  report  from  that  of  Allen 
lies  principally  in  the  evaluation  of  the  parameters  of  equations 
(10)  and  of  the  potential  theory.  It  was  well  known  to  both  Allen 

and  Hill  that  the  coefficient  CD  should  not  be  the  steady  state 

c 

drag  coefficient,  as  used  by  them,  but  should  be  related  to  the 

transient  effect  found  by  Schwabe  (Ref.  8).  The  results  of 

Schwabe's  experiments  are  reproduced  here  as  Fig,  2.  This 

curve  represents  the  drag  coefficient  for  a circular  cylinder 

moving  crosswise  in  a fluid  when  started  impulsively  from  rest. 

as  computed  from  the  observed  flows.  For  a cylinder  of 

radius  r,  , it  is  found  that  the  drag  coefficient  rises  rapidly  from 

V t 

zero  until  the  parameter  c is  about  4,  A very  slow  increase 

rl 

continues  until,  when  this  parameter  reaches  a value  of  9,  the 
drag  coefficient  is  2.07,  or  twice  the  steady-state  value  of  about 
1.0  for  the  Reynolds  number  of  the  experiment. 
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Referring  to  Fig.  1 we  »nay  consider  a plane  lamina  of  air 
of  thickness  dx  moving  along  the  body  axis  with  speed  U cos  a 
and  across  the  body  with  speed  U sin  a . At  the  nose  of  the  body, 
this  lamina  suddenly  begins  to  flow  about  a circular  cylinder  of 
variable  radius  r.  A potential  cross  flow  will  exist  at  the  nose, 
with  a gradual  development  of  viscous  cross  flow  along  the  body. 
Schwabe's  parameter  becomes 

Vc  U sin  ax  x 

rt_  r U cos  a r ‘an  a 

In  order  that  equations  (10)  may  be  integrated  readily,  it  is 

advisable  to  approximate  Schwabe's  results  by  a simple  function. 

This  function  was  first  chosen  to  be  a linear  increase  of  drag 

V_t  Vct 

coefficient  from  0 at  - 0 to  2.0  at  = 4,  and  to  remain 

r r 

V t 

constant  at  2.0  until  c = 9 . No  information  is  given  by 

r 

Schwabe  of  the  transition  back  to  the  steady  state  for  larger 
values  of  his  parameter.  It  may  be  seen  from  Fig.  3,  however, 
that  this  parameter  range  covers  most  body  lengths  and  angles 
of  Pttack  of  practical  interest. 
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The  linear  approximation  works  quite  well,  but  yields  un- 
wieldy expressions  for  angles  greater  than  that  for  which 
£ tan  a = 2,  since  the  factor  tan  a appears  as  a limit  of  integra- 
tion. The  experimental  curve  by  Schwabe  resembles  a portion  of 
an  arc  tangent  curve,  which  suggests  the  use  of  such  a function. 
The  evaluation  of  the  integrals  obtained  appears  to  yield  more 
complicated  formulas  than  the  present  problem  deserves.  The 
Taylor  expansion  of  an  arc  tangent  function 


-1 


y = tan  x = x 


3 + 5 + 


(12) 


suggests  that  the  same  result  may  be  obtained  by  integrating 
term  by  term,  except  that  the  factors  containing  tan  0 separate 
out  immediately  this  way.  The  limits  of  applicability  of  this 
theory,  due  to  a change  in  flow  conditions  at  large  angles,  renders 
the  use  of  terms  beyond  the  fifth  power  questionable,  so  that  a 
truncated  arc  tangent  of  the  form 


x3  i 


Cq  = A'  — tana+B'  — jtan3  a + C ' — ^-tan^  a + - - - (13) 


will  be  sufficient  for  practical  purposes. 
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The  cylinder  cross  flow  drag  function  may  be  integrated  for 
specific  body  shapes.  Since  the  integral  is  multiplied  by  sin  a, 
we  may  simplify  the  results  by  using 


sin  a tan  a = a- 


flZ 

15 


+ - - 


(14a) 


sin2  a 


tan3  a 


(14b) 


sin2  a tan*  a = a7  + - - - (14c) 

so  that,  to  terms  in  a7  , the  normal  force  and  moment  coefficient 
increments  will  be  of  the  form 

iCDc 

Acn  = -x— [DNa3  + EN05  + FN07)  (15a) 


2Cd 

AcM  = — A~[DMa3  + EMa?  + EMa?l  (l  5b) 


where  is  the  steady-state  value  of  the  cylinder  cross  flow 

drag  coefficient.  The  particular  functions  obtained  for  some 
typical  body  shapes  will  be  derived  in  Appendix  A.  It  has  been 


16 


TM  -998 


found  that  the  linear  term  in  Eq.  13,  leading  to  a cubic  in  normal 
force  coefficient  and  moment  coefficient,  is  sufficiently  accurate 
for  much  of  the  data  in  this  report. 

It  is  also  necessary  to  choose  an  appropriate  value  of  the 
parameter  T]  for  the  transient  phenomenon  being  considered.  The 
fact  that  H is  different  from  unity  for  the  steady  state  is  prin  • 
cipally  due  to  the  distortion  of  flow  near  the  ends  of  a cylindrical 
body.  A lamina  of  fluid  will  not  remain  in  a plane,  but  will  flow 
in  such  a way  that  the  total  flow  will  meet  less  resistance.  In  the 
case  of  the  transience velopment  being  considered  for  a blunt- 
based  body  of  revolution,  this  effect  will  not  have  time  to  develop 
at  the  base  before  the  flow  leaves  the  body.  A distortion  may 
develop  at  the  nose,  but  the  viscous  force  contribution  at  the  nose 
is  very  small,  and  the  effect  may  usually  be  neglected.  For  these 
reasons,  the  factor  rj  will  be  replaced  by  unity. 

If  a body  does  not  have  a blunt  base,  that  is  has  boattail,  the 
above  argument  will  fail,  since  considerable  distortion  of  the 
crods  flow  planes  will  occur  in  a region  where  the  viscous  effects 
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are  greatest  For  this  case,  the  value  of  T)  may  differ  signifi- 
cantly from  unity,  which  is  os'.e  reason  the  present  report  is 
restricted  to  blunt-based  bodies. 

An  interesting  feature  of  the  results  in  Eq.  15  is  that  the 
expressions  for  the  coefficients  become  odd  functions,  so  apply 
directly  to  both  positive  and  negative  angles.  Allen's  quadratic 
formulas  can  only  apply  to  negative  angles  of  attack  if  they  are 
considered  to  represent  absolute  values  of  the  coefficients. 

There  remains  the  problem  of  determining  the  appropriate 
value  of  (k  2 - ki  ),  the  factor  in  Munk'i  theory  representing  the 
difference  in  apparent  mass  for  longitudinal  and  transverse 
motion  of  the  body.  Munk  includes  a table  of  values  of  k2  - kt  , 
originally  computed  by  Lamb  (Ref.  9,  p.  155).  These  values 
were  computed  for  ellipsoids  of  revolution,  but  do  not  neces- 
sarily apply  to  cone  cylinder  or  ogive  cylinder  combinations  of 
the  same  fineness  ratio.  Munk  states  that  the  values  may  be 
used  for  airship  design  if  the  ellipsoid  has  the  same  length  and 
volume,  or  if  Volume^  has  the  same  value.  This  will  not  be 
as  good  an  approximation  for  a cone -cylinder , but  in  the  absence 
of  better  methods  will  serve  as  a method  for  correcting  what  is 
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already  a correction  factor.  A*  an  example  of  the  difference  this 
correction  makes  for  short  bodies,  a cone  cylinder  of  5 caliber 
total  length  and  2 caliber  head  length  is  found  to  be  equivalent  to 
an  ellipsoid  of  11  calibers  total  length.  The  value  of  k2  - k1  for 
5 calibers  is  0.836,  while  fcr  11  calibers  it  is  0.95.  Experimental 
evidence  favors  the  latter  v?lue. 

Boundary  Layer  Corrections 

Hill  has  pointed  out  in  Ref.  7 that  the  boundary  layer  may 
necessitate  two  corrections  to  the  theory  of  Allen,  one  in  the 
viscous  flow  term,  and  one  in  the  potential  term.  His  modified 
form  of  Allen's  theory  does  not  eliminate  the  fundamental  difficulty 
and  principal  source  of  error,  but  this  type  of  correction  may  not 
in  general  be  neglected  in  a precise  theory. 

Alien'll  viscous  cross  flow  term  assumes  that  the  cross  flow 
is  completely  independent  of  the  axial  flow  and  will  therefore 
depend  only  on  the  cross  flow  Reynolds  number  (Eq.  7a)  and  Mach 
Number  (Eq.  7b).  There  is  good  experimental  evidence,  how- 
ever (Ref.  6,  p.  431),  that  the  critical  Reynolds  number  for  the 
decrease  of  the  drag  coefficient  of  a cylinder  from  the  laminar 
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flow  value  to  the  turbulent  flow  value  changes  rapidly  with  the 

turbulence  level  present  in  the  stream.  It  seems  logical  that  if 

a cross  flow  in  this  range  of  Reynolds  numbers  encounters  a 

turbulent  axial  boundary  layer,  that  the  proper  drag  coefficient 

Cn  to  use  will  be  that  for  a turbulent  cross  flow.  The  coefficient 
uc 

appears  to  change  by  a factor  of  about  3,  and  the  steady  state 
values  of  Cpc  that  will  be  used  here  will  be  1 . 2 for  the  case  of 
completely  laminar  boundary  layer,  and  0.  35  for  (almost)  com- 
pletely turbulent  boundary  layer. 

The  effect  of  boundary  layer  on  the  potential  theory  term 
arises  from  the  fact  that  the  flow  is  displaced  as  though  the  body 
shape  were  changed.  In  particular,  the  effective  base  area  Sg 
and  volume  V will  be  greater.  The  important  boundary  layer 
parameter  will  be  the  displacement  thickness 


b 
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For  laminar  boundary  layer  thi«  thickness  may  be  found  approx 
imately  from 


b*  A 
* ' sTTT 

where 

A = 1.73  + 0 . 4 8M  2 


(17) 


(18) 


The  turbulent  boundary  layer  displacement  thickness,  as  shown 
by  Prandtl  from  the  seventh  power  law 


(19) 


is  given  by 


6 * 
x 


K 

0.2 

(Re) 


(20) 


where  K - 0.046  for  incompressible  flow  on  a flat  plate.  Flat 
plate  data  from  the  Daingerfield  supersonic  wind  tunnel  indicate 
an  approximate  value  of  K = 0.  08  for  a Mach  number  of  1. 56, 
whicn  is  of  interest  in  this  report.  Hill  has  computed  the 
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fractional  increase  in  base  area  from  Sg  to  Sg'  and  approximate 
formulas  for  the  fractional  increase  in  volume  from  V to  V', 
based  on  the  assumption  of  constant  diameter  and  fully  turbulent 
boundary  layer.  These  are,  for  the  laminar  case 


SB' 
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and  for  the  turbulent  case 
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Supersonic  Potential  Flow 

Allen  has  assumed  that  the  normal  force  and  pitching  moment 
may  be  found  by  the  simple  addition  of  the  potential  flow  solution 
and  the  viscous  cross  flow  solution.  The  logical  application  ol 
this  principle  will  use  the  best  available  potential  flow  solution. 
For  subsonic  flow,  the  theory  of  Munk,  corrected  for  boundary 
layer,  gives  excellent  results.  For  supersonic  flow,  however, 
Munk's  theory  does  not  apply,  except  possibly  in  the  limiting  case 
ot  bodies  with  extremely  long,  slender  noses. 

A much  better  approximation  to  the  potential  flow  solution 
tor  supersonic  flows  is  obtained  by  the  second-order  theory  of 
V in  Dyke  (Ref.  10).  This  theory  has  been  applied  recently  to  the 
calculation  of  pressure  distributions,  normal  force  coefficients, 
and  pitching  moment  coefficients  for  a series  of  tangent -ogive  - 
Cylinder  combinations  (Ref.  11).  These  must  necessarily  be 
computed  for  a finite  number  of  discrete  supersonic  Mach 
numbers,  so  a method  of  interpolation  is  desirable  In  Fig.  4 
and  5.  the  values  of  dCjy/da  and  dCj^j/da  for  a = 0 degrees 


as  computed  m Ref.  11  are  plotted  as  a function  of  j3  - si  _ j 
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MACH  PARAMETER,  p = M2  - 1 

FIG.  4.  INITIAL  SLOPE  OF  COEFFICIENT  OF 
PITCHING  MOMENT  (ABOUT  NOSE)  FROM 
VAN  DYKE  1HEORY,  j£  = 14- 
14 
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3.5  CALIBER  HEAD  LENGTH 


2.5  CALIBER  HEAD  LENGTH 


■■■■■■■■■■■■ 

■■■■■■■■■■ 


MACH  PARAMETER,  P ^ \2  - i 

. INITIAL  SLUPE  OF  COEFFICIENT  OF  NORMAL 
FORCE  FROM  VAN  DYKE  THEORY,  JL  =14. 
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with  ogival  head  length  as  a parameter.  These  values  were 
computed  tor  a total  body  length  of  14  calibers,  and  do  not  change 
appreciably  with  increase  in  body  length  beyond  9 calibers.  It 
may  be  seen  that  interpolation  may  be  made  along  smooth  curves 
for  the  moment  coefficient  slope,  while  the  interpolation  is  linear 
for  the  normal  force  coefficient  slope.  No  particular  significance 
has  been  attached  to  this  linear  relation,  except  for  a convenient 
method  of  interpolation. 

EXPERIMENTAL  VERIFICATION 
The  theory  described  above  is  compared  with  Allen's  method 
and  with  experimental  data  in  Appendix  B,  Fig.  6-19.  The  extra- 
polation of  Munk's  linear  potential  theory  in  the  form  Cj^  = Z a 
2 V 

and  Cjy(n  = — — Ct  is  included  as  a reference  line  for  demonstration 

A 

of  nonlinearity.  The  exact  body  shape  was  not  used  in  computing 
theoretical  results.  For  the  tangent  ogives,  the  circular  arc  was 
approximated  by  a parabolic  arc,  as  explained  in  Appendix  A.  The 
secant  ogives  for  the  short  bodies  were  approximated  by  rones. 

The  contribution  to  the  viscous  cross  force  from  the  head  is  so 
small  that  the  error  involved  in  these  approximations  is  negligible. 
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Short  Bodies  in  Subsonic  Flow,  With  Turbulent  Bounfary  Layer 

The  moment  coefficient  for  the  5 caliber  body  (Fig.  6)  shows 
very  little  nonlinearity  up  to  20  degrees,  but  in  the  range  from  20 
to  30  degrees  a slight  curvature  is  shown.  The  agreement  with 
experimental  data  is  everywhere  excellent.  Allen's  method 
grossly  overestimates  the  nonlinear  effects.  The  normal  force 
coefficient  vFig.'  7)  shows  a similar  agreement  between  the 
present  theory  and  experiment,  with  overestimation  by  Allen's 
method.  No  definite  conclusions  may  be  drawn  from  the  center  of 
pressure,  except  that  both  nonlinear  theories  predict  the  rearward 
shift  with  large  angles  of  attack. 

The  moment  coefficient  (Fig.  8)  and  the  normal  force  co 
efficient  (Fig.  9)  for  the  7 caliber  body  again  show  an  excellent 
agreement  between  the  present  theory  and  experiment,  Allen's 
method  again  being  high.  The  center  of  pressure  shows  a larger 
rearward  shift  with  angle,  with  slightly  better  agreement  for  the 
present  theory  than  for  Allen's  method. 

The  correlation  of  the  present  theory  with  experiment  for  the 
9 caliber  length  is  also  remarkably  good  (Fig.  10,  1 1 ) . The 
results  tor  center  of  pressure  are  more  distinct  in  this  case. 
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Comparing  the  low-angle  results  here  with  those  in  Fig.  7 and 
Fig  9,  it  seems  likely  that  the  boundary  layer  corrections  are 
inadequate.  This  deviation  is  not  serious,  however,  considering 
the  ser.si  tivity  of  center  of  pressure  to  small  errors  in  forces 
and  moments . 

It  is  noteworthy  that  Allen's  method,  based  on  subsonic 
theory  and  data,  shows  a greater  deviation  from  experiment  at 
these  low  speeds  than  for  low  supersonic  speeds.  It  may  be 
argued  that  it  is  not  a fair  comparison  to  correct  for  the  apparent 
mast  factor,  boundary  layer  thickness,  and  for  the  fact  that 
sin  L a differs  from  2a  in  presenting  the  results  of  the  present 
theory,  but  not  in  Allen's  formula.  Such  corrections  in  Allen's 
formula,  however,  would  account  for  only  a small  part  of  the 
error  shown,  and  the  formula  is  used  just  as  used  for  large  angles 
by  Allen  and  Perkins  (Ref.  4). 

Turbulent  Boundary  Layer  at  Supersonic  Speed 

Data  are  presented  for  three  different  bodies  of  revolution  at 
a Mach  number  M = 1.  56  in  Fig.  12  17.  The  Reynolds  numbers 
for  these  tests  were  such  that  the  boundary  layer  was  probably 
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turbulent  for  must  of  the  body  length,  except  for  the  10  caliber 
model  Schlieren  pictures  have  indicated  a boundary  layer  transi- 
tion about  5 calibers  from  the  nose. 

The  agreement  with  experiment  is  better  on  the  average  for 
the  present  theory  than  for  Allen's  method.  The  apparent  agree 
rnrnt  between  his  method  and  the  data  is  probably  fortuitous,  being 
the  result  of  two  or  more  errors  tending  to  cancel.  This  is  espe- 
cially true  for  the  longer  bodies,  where  the  initial  slope  for  the 
experimental  data  is  considerably  greater  than  for  Allen's  theory, 
which  has  the  same  initial  slope  as  the  linear  potential  theory. 

The  agreement  between  the  present  theory  and  experiment  is 
good  for  center  of  pressure,  where  deviations  of  the  order  of  half 
a caliber  may  be  within  experimental  error. 

Laminar  Boundary  Layer  at  Supersonic  Speed 

One  body  shape  has  been  selected  for  which  the  boundary 
layer  probably  remains  laminar  at  the  test  Mach  number  of  2.87 
and  Reynolds  number  of  3.72  million.  It  was  also  chosen  because 
Van  Dyke's  second  order  theory  applies  to  it  at  this  high  Mach 
number  For  this  body,  of  length  14  calibers  and  head  length 


29 


TM  -9°8 


1.5  calibers  the  corrections  for  displacement  ‘Sickness  of  the 
laminar  boundary  layer  are  small  compared  to  the  turbulent 
thickness  at  Mach  1.56,  but  are  included  for  completeness.  The 
value  Cn  - 1 L is  used  in  contrast  to  Cn  = 0.35  for  all  the 
previous  curves. 

The  failure  of  Alien's  formula  is  quite  apparent  in  Fig.  18 
and  19  For  the  true  laminar  case  shown,  his  combination  of  sub- 
sonic potential  theory  and  quadratic  viscous  term  falls  far  short 
of  the  experimental  data,  even  for  angles  up  to  8 degrees.  A 
comparison  of  Fig.  18  and  19  with  Fig.  16  and  17  shows  the  striking 
difference  in  nonlinear  effects  between  laminar  and  turbulent 
h undar y layers . 

PRACTICAL  ENGINEERING  APPLICATION 

It  has  been  shown  that  the  present  theory  agrees  well  with 
experimental  data  over  a wide  range  of  Mach  numbers,  Reynolds 
numbers  and  body  shapes.  For  practical  use,  it  is  desirable  to 
be  able  to  simplify  the  equations  whenever  possible.  The  complete 
equations  may  be  given  as 
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dCN  dcM 

The  and  — yu—  are  given  by  Munk's  theory  at  subsonic 

da  du 

speeds,  after  correcting  for  the  apparent  mass  factor  and  the 
boundary  layer  thickness,  and  for  supersonic  speeds  are  found 
from  Van  Dyke's  second  order  theory,  with  boundary  layer  cor- 
rections added.  The  D,  E,  F coefficients  are  computed  as  shown 
in  Appendix  A,  ami  is  1.2  or  0.35  as  the  axial  boundary  layer 

is  laminar  or  turbulent,  respectively. 

Figure  20  shows  a comparison  of  the  normal  force  coefficient 
Cj^j,  as  computed  by  the  complete  equation  with  various  approxi- 
mations, for  the  9 caliber  length  body  at  subsonic  speeds.  It  is 
found  that  the  combined  corrections  for  boundary  layer  thickness 
and  apparent  mass  factor  may  be  omitted  for  this  case  with  little 
loss  in  accuracy.  These  corrections  tend  to  cancel,  and  are  of 
the  same  o^der  of  magnitude  for  body  lengths  of  5 to  7 calibers. 
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For  very  nhort  bodies  the  (k2  - 1c,  ) factor  will  dominate  and  for 
very  long  bodies,  with  k2  - k,=  1 the  boundary  layer  correction 
will  become  important. 

The  use  of  2 sin  a cos  a cos  y ? sin  2a  cos  y **  2a  is  shown 
to  be  valid  up  to  10  degrees,  and  has  been  used  this  way  for  all 
figures  for  which  the  angle  of  attack  a was  less  than  8 degrees. 
Approximations  obtained  by  dropping  the  seventh  power  term,  or 
both  fift'i  and  seventh  powers,  show  that  the  cubic  equation  is 
quite  accurate  for  2 tana  3,  the  quintic  equation  holds  to 
iJl  tan  a = 6,  and  the  seventh  power  is  good  to  2 J^tan  a = 9,  the 
limit  of  application  of  this  theory  (see  Fig.  3). 

The  case  of  turbulent  axial  boundary  layer  in  supersonic  flow 
is  shown  in  Fig.  21.  Normal  force  coefficients  are  plotted  for  a 
14  caliber  length  body  with  2.8  caliber  nose  length  for  the  complete 
equation  (but  with  sin  2 a = 2a)  and  with  approximate  formulas. 

The  neglect  of  boundary  layer  thickness  begins  to  be  significant, 
if  an  accurate  estimate  of  normal  force  is  desired,  but  may  be 
sufficiently  good  for  rough  estimation  purposes.  The  curves  are 
approximated  very  well  by  a cubic  equation,  in  agreement  with 
the  criterion  of  the  last  paragraph.  The  errors  due  to  boundary 
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layer  thickness  will  increase  with  Mach  number  and  body  length, 
and  the  errors  due  to  shortening  the  polynomial  expression  in- 
crease with  body  length,  for  a given  angle  of  attack. 

Comparisons  are  also  made  for  the  case  of  laminar  boundary 
layer  (Fig.  22)  on  a 14  caliber  body.  The  results  of  shortening 
the  polynomial  are  similar  to  those  found  in  the  last  paragraph. 
For  the  laminar  boundary  layer,  the  effects  of  displacement 
thickness  are  small  in  comparison  with  similar  effects  of  a 
turbulent  boundary  layer,  and  may  usually  be  neglected. 

It  is  found  that  for  a special  case  of  Bubjonic  flow,  with 
moderate  length  bodies  and  a i 10  degrees,  the  formulas 


CD  J 

CN  = 2 a * 0.49— - — a3 


(26a) 


Cmb  = ^“+  o.« 


cDci' 

3A 


(26b) 


where  Dj^  = 0.49  ^ , Dj^  = 0.49^  have  been  substituted, 

give  good  approximations  to  the  normal  force  and  pitching 
moment  coefficients.  These  formulas  will  not  always  be  suffi- 
cient for  large  angles,  very  long  bodies,  for  supersonic  flows. 
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or  for  the  rase  of  very  thick  turbulent  boundary  layer,  but  the 

corrections  for  these  cases  are  usually  simple  to  apply.  The 

value  of  must  be  chosen  as  1.2  or  0.35  as  the  axial  boundary 

layer  is  laminar  or  turbulent,  respectively. 

In  the  case  of  supersonic  flow,  there  are  some  body  shapes 

and  Mach  numbers  for  which  Van  Dyke's  theory  does  not  apply 

This  theory  is  still  useful  if  experimental  data  at  small  angles 

dCN  dCj^j 

of  attack  are  available.  Slopes  and  — . at  a = 0 may 

da  da 

be  measured  and  used  in  Eq  25  in  place  of  the  slopes  from  Van 
Dyke's  theory.  This  technique  will  be  valuable  in  extrapolating 
wind  tunnel  data  for  small  models  to  free-fltght  conditions  with 
large  models  Care  should  be  taken  to  note  the  boundary  layer 
conditions  of  the  wind  tunnel  tests  and  free  flight  applications. 

CONCLUSION 

A method  has  been  developed  for  the  estimation  of  normal 
force  coefficients  and  ptching  moment  coefficients  for  blunt- 
based  bodies  of  revolution  at  large  yaws  The  equations  reduce 
to 
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CD  fi-2 

CN  = 2at  0.49  — a3 
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mb 
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for  subsonic  speeds  if  a £ 10  degrees,  2j^tanC(  £ 3 , as  long  as  the 
boundary  layer  thickness  does  not  increase  the  base  area  appre 
ciably.  Van  Dyke's  theory  should  be  used  for  the  linear  term  in 
supersonic  flow,  and  a higher  order  polynomial  is  needed  if 
zltand  2 3 . 

The  method  is  also  useful  for  extrapolation  of  small-yaw 
data  to  large  yaws  and  to  different  Reynolds  numbers.  As  far  as 
is  known,  this  is  the  first  adequate  theory  that  has  been  developed 
for  this  purpose. 

The  results  presented  here  have  only  been  applied  in  the  Mach 
number  range  from  M = 0 to  M = 2.87  and  for  a limited  range  c.'. 
Reynolds  numbers.  Further  investigation  will  be  made,  especially 
in  the  high  Mach  number  range.  The  results  of  a more  complete 
investigation  will  be  described  in  a later  report. 
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Appendix  A 

CALCULATION  OF  NONLINEAR  TERMS 

The  experimental  measurement  by  Schwabe  of  the  development 
of  crons  drag  coefficient  for  a circular  cylinder  (Fig.  Z)  resembles 
an  arc  tangent  function,  except  at  small  values  of  the  argument, 
where  contributions  to  normal  force  in  the  present  theory  will 
usually  be  quite  small  in  comparison  with  the  linear  term.  The  arc 
tangent  function  may  be  integrated  directly,  but  results  in  compli- 
cated formulas.  A Taylor  expansion  of  the  arc  tangent,  which  is 
an  infinite  series  containing  only  odd  powers,  converges  quite 
slowly.  Such  a series  may  be  truncated,  however,  so  that  it  rep- 
resents the  function  over  a limited  range  of  the  argument  with  the 
use  of  a few  of  the  leading  terms,  if  the  coefficients  of  these  terms 
are  properly  adjusted.  The  representation  of  Schwabe's  data  was 
chosen  to  be  of  the  form 

Cn  = A'- tana  + B 1 — r-  tan3Cl  + C'^r-tan5a  (lA) 

uc  r r3  r? 
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By  a least  squares  procedure,  A',  B1,  and  C!  were  evaluated 
giving 

CDc  = 0.49^-tana  - 0.  0056  l^.tan3a  + 0.  00003  2^  tan5a  (2A) 


This  function  is  compared  with  Schwabe's  data  in  Fig  2. 

For  computation  of  viscous  cross  flow  effects,  one  must 
evaluate  the  integrals  in  the  expressions 


AC 


N 1 r ,in2a / 


r CDc  dx 


Ac 


M 


= j.  .in ‘a/ 

•'O 


r CD  (xm  * x)dx 


(3A) 


(4A) 


where  C^  is  given  by  ( 1 A ) multiplied  by  the  steady-state  value, 
since  the  steady  state  value  in  Schwabe's  experiment  was  about 
unity.  This  procedure  assumes  that  the  shape  of  the  curve  found 
by  Schwabe  always  applies,  and  only  the  absolute  magnitudes 
change  with  Reynolds  number.  These  integrations  have  been 
carried  out  for  bodies  with  conical  and  parabolic  ogival  heads  on 
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cylindrical  bodies,  where  the  body  length  is  JC  calibers  and  the 
head  length  is  h calibers.  The  radius  r for  the  generating 
curves  of  the  heads  are  given  by 


r 


(5A) 


calibers  for  the  cone  and 


(6A) 


for  the  parabola.  The  parabola  is  a close  approximation  to  a 
circular  tangent  ogive,  and  is  much  simpler  to  integrate  than  the 
exact  function. 

It  is  sufficient  to  perform  the  integrations  for  the  case  of  a 
steady  state  Cr\  of  unity,  and  apply  a correction  factor  outside 
the  integral.  The  integral  for  the  normal  force  coefficient  for  a 
cone  -cylinder  is 
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dx  = A'  — - tan 


a + B'UU 


+ hu  Jtan^a 


+ 16C[ 


i6 

T 


u6 

-j— ]tan^a 


(?A) 


39 


TM  -998 


and  for  a parabolic -ogive -cylinder 


+ + + log  Zjtan^a 


+ 


16C'[ 


~T 


(8A) 


For  the  moment  coefficients  about  a point  xm  calibers  from  the 
nose  of  he  body,  the  integral  for  the  cone -cylinder  is 


/*£  2x 

/r  CDc(xm  - x)dx  = - ^]tand  + B '[(jgu  f hU 


lm 


/4i5  8h5.,  , , .,0*  h6 

' 5 ~i3“)]tan  a+  l6C'[^g—  + — )x 

<4- + Tr,),,n5a 


6 h6 

3~^m 


(9A) 


and  for  the  ogive  -cylinder 


Dc(xm  - x)dx  = A'[~  m .^]tana+  B'[(je4+  ~~  + log  2)xm 

/4  £5  36h5  , . . p6 

1“  ~T~  + l6h  io8  h)]tan3a  + 16C’[(^-  +|j)xm 

' * Tz  ^tan5<X  ( 1 0 A ) 
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These  expressions  simplify  for  the  case  xm  - Ji  , for  the 
moment  coefficient  to 


+ 16C'[^  + -^"Itan3 


(11A) 


for  the  cone,  and 


for  the  ogive. 

The  powers  of  tan  Q will  combine  with  sin2a  to  give,  to 
seventh  degree  terms. 


sin2a  tan  a = a3  + Sil.  + - 

15 


sin*  a tan3a  - ^ — + - - - 


2a1 


(13A) 

(HA) 
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sin2  a tan1’ a = a7  + 


(15A) 


These  may  be  substituted  in  Eq . 3A,  4A,  7A,  and  11A  to  give, 
for  cone-cylinders, 

2CDc  A'i2 

iCN  = -A-£l:iT-“3+  h*  )a* 

+ (16C  |4"  + t]  + ^U+h4}  + ^£2)a7]  (16A) 


Ac 


2C, 


Mg  " a 


(A£a3,B,,iltihV!*V 


*’**'-£} 


♦&iV  j 


(17A) 


where  Cn  is  now  the  steady-state  value  of  the  cylinder  cross 
c 

flow  drag  coefficient. 

A similar  substitution  for  ogive -cylinde r s yields  from  Eq. 
3 A,  4 A,  8A,  and  12A 
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(ISA) 


Acmb  = •2^r£t^“  °3  + B'^  + ^r~  + * 108  2 -3t- 

- I6h  log  h)a5  + (16C‘  ^ + 4—  + 5t|  +3"  fer 

+ ^ 1°8  2 - - I6h  log  hj  + ^p)*7]  (19A) 

Many  of  the  terms  in  16A  through  19A  may  be  neglected, 
especially  for  h<  <J?. 

As  an  example  of  the  use  of  the  above  equations,  the  viscous 
increments  for  a 7 caliber  cone -cylinder , j£  = 7 and  h = 2,  become, 
with  C = 0.35 

Acn  = pp  [l2a3  - 13.  5a5  + I.  22a7]  (20A) 


4 3 
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y ft 

AcMB  * ^-[28a3  - 19.  35a5  -l.55a7]  (21A) 

The  ogive -cylinder  with  Jl  - 14,  h = 3.8,  gives 

A CN  = ^[48a  3 - 216a  5 + 460a7]  (22A) 

AcMfi  - i^i[224a3  - 591a5  + 824a7]  (23A) 

9.  6 

The  front  multiplying  factor  must  be  changed  to  -75^—  for  a 
laminar  boundary  layer. 
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COMPARISON  OF  THEORIES  WITH 
EXPERIMENTAL  DATA 
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a,  ANGLE  OF  YAW,  DEGREES 


FIG.  6.  PITCHING  HOMENT  COEFFICIENT  FOR  A 
SECANT-OGIVE-CYLINDER  AT  U = 0.26. 

5 CALIBER  LENGTH 
2 CALIBER  HEAD  LENGTH 

Re  - 4.2  x 105 6 


NORMAL  FORCE  COEFFICIERT 
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FIG.  7.  NORMAL  FORCE  COEFFICIENT  AND  CENTER  OF 
PRESSURE  FOR  A TEC ANT- OGIVE-CYLINDER 
AT  M = 0.26. 


5 CALIBER  LENGTH 
2 CALIBER  HEAD  LENGTH 

Re  = 4.2  x 105  6 
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— THEORY  OP'  THIS  REFORT 

— ALLEN'S  FORMULA  (EQ.  5) 

— UNEAR  POTENTIAL  THEORY 
O EXPERIMENTAL  DATA 


VI 

- 

0 10  20 
a,  ANGLE  OF  YAW,  DEGREES 

FIG.  9.  NORMAL  FORCE  COEFFICIENT  AND  CENTER 
CF  PRESSURE  FOR  A SECANT-OGIVE-CYLINDER 
AT  M = 0.26 

7 CALIBER  LENGTH 
2 CAUBER  HEAD  LENGTH 
Re  = 5.9  x 106 


CENTER  CF  PRESSURE,  CALIBERS 


MO.  'JPT  COEFFICIENT  ABOUT  BASE 
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— THEORY  OF  THll  REPORT 

— ALLEN'S  FORMUB.  \ (EQ.  5) 

— LINEAR  POTENTU  L,  THEORY 
O EXPERIMENTAL  DAV'A 
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Z2 

' yS 

^4 

\ 

. 

1 

a angle  of  yaw,  degrees  \ 

FIG.  11.  NOR.  \L  FORCE  COEFFICIENT  AND  CEN1\H 
OF  FRESSUR1  FOR  A SEC ANT- OGIVE-CYLINDER  1 
AT  M = 0.26.  \» 

\ 

9 i \ LIBER  LENGTH 

2 C: LIBER  HEAD  LENGTH  |j 

Re  - 7.6  x TO6  - 


SI 


CENTER  CF  FRETS  LIRE,  CALIBERS  FROM  BATE 
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FIG.  12.  PITCHING  MOMENT  COEFFICIENT  FOR  A 
TANGENT-OGIVE-CYLINDER  AT  M = 1.56. 


10  CALIBER  LENGTH 
2.5  CALIBER  HEAD  LENGTH 
Re  = 4 x 10^ 
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FIG.  13.  NORMAL  FORCE  COEFFICIENT  AND  CENTER 
OF  PRESSURE  FOR  A TANGENT- OGIVE-CYLINDER 
AT  U = 1.56. 


10  CALIBER  LENGTH 
2.5  CALIBER  HEAD  LENGTH 
Re  = 4 x 10 6 
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a,  ANGLE  OF  ATTACK,  DECREES 


FIG.  H.  PITCHING  MOMENT  COEFFICIENT  F 
TANQBNT-OQIVE-CTUNDER  AT  M = 1.56. 


14  CALIBER  LENGTH 
2.5  CALIBER  HEAD  LENGTH 
Re  = 5.6  x 106 


r 
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CENTER  CF  PRESSURE,  CALIBERS  FRO”  BASE 
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- THEORY  OF  THIS  REPORT 

- ALLEN'S  FORMULA  (EQ.  5) 

- LINEAR  POTENTIAL  THEORY 
O EXPERIMENTAL  DATA 


mm 


a,  ANGLE  OF  ATTACK,  DEGREES 

FIG.  17.  NORMAL  FORCE  COEFFICIENT  AND  CENTER 
OF'  PRESSURE  FOR  A TANGENT- OGIVE-CYLINDER 
AT  M = 1.56. 

14  CALIBERS  LENGTH 
3.5  CALIBERS  HEAD  LENGTH 


5.6  x 10 
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a,  ANGLE  OF  ATTACK,  DEGREES 


FIG.  IS.  PITCHING  MOMENT  COEFFICIENT  FOR  A 
TANGENT- OGIVE-CYLINDER  AT  M - 2.87  , 

U CALIBER  LENGTH 
3.5  CALIBER  HEAD  LENGTH 
Re  = 3.7  x 106 


NORMAL  FORCE  COEFFICIENT 
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FIG.  19.  NORMAL  FORCE  COEFFICIENT  AND  CENTER 
QF  PRESSURE  FOR  A TANGENT- OGIVE-CYLINDER 
AT  M » 2.37. 

H CALIBER  LENGTH 
3.5  CALIBER  HEAD  LENGTH 

Re  = 3.7  x 106 
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a , ANGLE  OF  YAW,  DEGREES 


FIG.  20.  NORMAL  FORCE  COEFFICIENT  FOR  A 
NINE  CALIBER  BODY,  USING  DIFFERENT 
APPROXIMATIONS „ 
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NOMENCLATURE 


A 

Ap 

A',  B',  C' 


Dc 


^ L 
CM 

CN 

D,  E,  F 

Dc 

K 

L 

M 

Mc 

Mp 

N 


Reference  area  for  defining  coefficients 
Planform  area  of  body 
Polynomial  coefficients  (Eq.  13) 


D, 


Drag  coefficient  of  cylinder  in  a crossflow  —-7 

PVA 

2 AP 


Lift  coefficient 


qA 


Pitching  moment  coefficient 


Mp 

qAd 


Pitching  moment  coefficient  about  base 

N 


Normal  force  coefficient 


qA 


Polynomial  coefficients  (Eq.  15) 

Crossflow  drag  on  a cylinder 

Parameter  for  turbulent  boundary  layer  (Eq.  20) 
Lift  force 
Mach  number 

Mach  number  of  crossflow 
Pitching  moment 
Normal  force 


Re 


Reynolds  number 
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Rec 


SB' 


u 

V 


V' 


Reynolds  number  of  crossflow 

Gr  js 8 - sectional  area  of  body  at  station  x 

Cross-sectional  area  at  base 

Base  area,  corrected  for  boundary  layer 

Free  stream  velocity 

Volume  of  body 

Volume,  corrected  for  boundary  layer 
Velocity  of  crossflow 
Speed  of  sound 

Diameter  of  body  ( = 1 caliber) 

Head  length  in  calibers 

Apparent  mass  correction  from  Munk's  theor 

Fineness  ratio  (body  length  in  calibers) 

n U* 

Dynamic  pressure  p — 

Body  radius  at  station  x 

Time  in  seconds 

Local  axial  velocity  in  boundary  layer 
Axial  coordinate  in  calibers 
Axial  position  of  body  centroid 
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Center  of  momenta 


Angle  of  attack 
^ M 2 - 1 

Boundary  layer  thickness 
Boundary  layer  displacement  thickness 
Paramek  • for  laminar  boundary  layer  (Eq.  17,  18) 
Drag  correction  factor  for  finite  length  cylinder 


Kinematic  viscosity  of  air 
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